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Abstract The leading hadronic contributions to the 
anomalous magnetic moments of the electron and the 
T-lepton are determined by a four-flavour lattice QCD 
computation with twisted mass fermions. The results 
presented are based on the quark-connected contribu¬ 
tion to the hadronic vacuum polarization function. The 
continuum limit is taken and systematic uncertainties 
are quantified. Full agreement with results obtained by 
phenomenological analyses is found. 
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1 Introduction 

The standard model of particle physics (SM) contains 
three charged leptons I, mainly differing in mass, the 
electron, the muon, and the r-lepton with rrie : : 

itIt ~ 1 : 207 : 3477 [I]. Their magnetic moments, 
in particular their so-called anomalous magnetic mo¬ 
ments, ai = {g — 2 ) 1 / 2 , control their behaviour in an 
external magnetic field. 

Being the lepton with the smallest mass, the elec¬ 
tron is stable. This leads to the electron magnetic mo¬ 
ment being one of the most precisely determined quan¬ 
tities in nature. The difference between the latest ex¬ 
perimental [5] and SM values [?1E] is of O (10 or 
approximately 1.3 standard deviations, c.f. and ref¬ 
erences therein, 
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= 115 965 218 07.3(2.8) x 10^^ 
a™ = 115 965 21817.8(7.6) x 10"^^ 
af’^P - a™ = -10.5 (8.1) x lO'^^ . 

This constitutes one of the cornerstone results for 
quantum field theories to be recognised as the correct 
mechanism for describing particle interactions. The very 
good agreement of the electron magnetic moment be¬ 
tween experiment and SM calculations is not matched 
by the muon anomalous magnetic moment. In fact, here 
a two to four sigma discrepancy is observed, see e.g. [6]. 
One reason for the observed discrepancy could be that 
the magnetic moment of the muon receives larger non - 
perturbative contributions than the one of the electron. 
On the other hand, it is supposed to be also more sensi¬ 
tive to beyond the SM physics, since for a large class of 
theories new physics contributions are expected to be 
proportional to the squared lepton mass. Thus it is a 
prime candidate for detecting physics beyond the SM. 
Due to the large mass of the r-lepton, it would be the 
optimal lepton for finding new physics. However, be¬ 
cause its lifetime is very short {O (lO“^^)s) there cur¬ 
rently only exist bounds on its anomalous magnetic mo¬ 
ment from indirect measurements [7]. 

The QED [^[3] and the electroweak contributions li 
m to the lepton anomalous magnetic moments have 
been computed in perturbation theory to impressive 
five and two loops, respectively. The main uncertain¬ 
ties remaining in the theoretical determinations of the 
anomalous magnetic moments originate thus from the 
leading-order (LO) hadronic contributions. Since they 
are particularly sensitive to those virtual photon mo¬ 
menta that are of O these contributions are in¬ 

herently non-perturbative and not accessible to per¬ 
turbation theory. In order to have a prediction of the 
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anomalous magnetic moments from the SM alone, a 
non-perturbative method needs to be employed and the 
only such approach we presently know, which eventu¬ 
ally allows us to control all systematic uncertainties, is 
lattice QCD (LQCD) which we use here. 

As highlighted in [^, the uncertainty in the com¬ 
parison between the experimental and the SM value for 
the electron anomalous magnetic moment is currently 
dominated by the experimental uncertainty of its deter¬ 
mination and the value for aqED from atomic physics 
experiments with rubidium atoms which both are to be 
reduced in the future. Recently, the Harvard group has 
announced to be working on a more accurate determi¬ 
nation of the electron as well as the positron {g — 2) [TT] . 
According to Ref. |S], uncertainties in the sub-10“^^ re¬ 
gion might be expected which would clearly provide the 
opportunity to also detect new physics contributions in 
the anomalous magnetic moment of the electron and 
thus to cross-check the muon discrepancy. In this situa¬ 
tion it will again be of utmost importance to know the 
hadronic contributions as precisely as possible. 

Furthermore, even for the r-lepton. Ref. m lists 
several proposals for the first actual measurement of its 
anomalous magnetic moment, e.g. m- A hrst success¬ 
ful measurement in this direction has been reported in 
M- As we will show in the following, compared to the 
case of the muon it will be much easier to obtain a value 
for the LO QCD contribution to Ur from LQCD with 
the required precision to detect new physics and it will 
probably not take very long before the QCD contribu¬ 
tion entering the official SM result will be provided by 
LQCD. 

As mentioned before, the hadronic LO contributions 
to the anomalous magnetic moments of the three SM 
leptons, strongly depend on the values of their 

masses. Since the magnitude of the lepton masses spans 
four orders of magnitude, the corresponding contribu¬ 
tions to the anomalous magnetic moments differ sub¬ 
stantially and probe very different energy regions, see 
also the discussion of Fig. in Sect. 

In this article, we present the results of our four- 
flavour computations of the quark-connected, LO hadronic 
vacuum polarisation contributions to the electron and 
T-lepton anomalous magnetic moments obtained from 
the (maximally) twisted mass formulation of LQCD. 
The muon case has already been covered in |15j . The 
important feature of the present calculation is that we 
adopt exactly the same strategy as for the muon m in¬ 
cluding the same chiral and continuum extrapolations. 
Thus, the results presented here are not only interest¬ 
ing in themselves, but also serve as an important cross¬ 
check for our treatment of the hadronic vacuum polar¬ 
isation function. The consistent picture we obtain for 


all three standard model leptons then reassures the va¬ 
lidity of the analysis approach. 

Additionally to the systematic uncertainties investi¬ 
gated in our previous paper, we quantify the light-quark 
disconnected contributions on one of our Nf = 2-l-H-l 
ensembles in order to arrive at a rough estimate of their 
systematic effect on our estimates for the hadronic LO 
anomalous magnetic moments. A full quantitative con¬ 
straint on the quark-disconnected contribution, how¬ 
ever, is beyond the scope of this work. 

Another very important feature is that incorporat¬ 
ing the complete first two generations of quarks en¬ 
ables us to directly and unambiguously compare our 
results with the values obtained from phenomenolog¬ 
ical analyses relying on experimental data and a dis¬ 
persion relation. We note that the contributions from 
third-generation quarks can be neglected, since they are 
smaller than the current theoretical accuracy, as can be 
inferred e.g. from the data tables of Ref. m- Recently, 
the bottom quark contribution to has been explic¬ 
itly computed on the lattice m confirming it to be one 
order of magnitude smaller than the current uncertainty 
of the phenomenological determinations of aJC^P. 

Additionally to the Nf = 2 -|- 1 -t- 1 flavour ensem¬ 
bles [T^fT^ at unphysically large pion masses studied 
in we computed the dominant light quark contribu¬ 
tions to the anomalous magnetic moments on a Nf = 2 
flavour ensemble directly at the physical point Ennu. 
This allows us to test the chiral extrapolations per¬ 
formed when using the reparametrisation introduced 
in mna. Since we currently only have one ensemble 
at the physical point with one lattice spacing and we 
neglect the small influence of the strange and charm 
sea quarks when comparing the results from Nf = 2 
and Nf = 2 -|- 1 -I- 1 ensembles, a final conclusion is pre¬ 
cluded at this point. The significance of this comparison 
is based on the empirically observed weak dependence 
on the lattice spacing of the light quark contribution as 
well as the marginal sea quark effects from strange and 
charm on the latter. 

The next section comprises a short repetition of the 
most important equations needed to follow the discus¬ 
sion of the results for the LO hadronic vacuum polarisa¬ 
tion contributions to the anomalous magnetic moments 
of the electron in Sect. and the r-lepton in Sect.|^ In 
Sect. [5] we summarise our results and draw our conclu¬ 
sions. 


2 Computation of 

The LO hadronic contribution to the lepton anomalous 
magnetic moments, can be directly computed in 
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Euclidean space-time according to (SUES] 

hvp 2rdQ^ m 

where a is the fine structure constant, the Euclidean 
momentum, mi the lepton mass, and n^{Q^) the renor¬ 
malised hadronic vacuum polarisation function, 

n^iQ^) = n{Q^) - iT(o). 

It is obtained from the hadronic vacuum polarisation 
tensor 

n^AQ) = J d‘^xe^^<^-yHj^{x)My)) 

= (Qi^Qu - n{Q ^), (2) 

which is transverse because of the conservation of the 
electromagnetic current 

2 1 - 

= i^u{x)-f^u{x) - - d{x)-f^d{x) 

2 1 

+ 3 c(a:)7^c(x) - - s{x)-f^s{x) . (3) 


Refs. [22ll2^ proposed a modified definition of a new 
family of observables 




mf ) 




( 5 ) 


H denotes some hadronic scale determined on the lat¬ 
tice at unphysically high pion masses, which fulfills the 
constraint, that H (mps) I^phys as mps —>■ TO,r- 
For each choice of H we thus obtain a correspondingly 
modified, mps dependent lepton mass on the lattice 
«r/-(mps) = mi ■ H (tops) /I?phys- 

Our choice for the hadronic scale H is the mass my 
of the lowest-lying state in the light vector meson chan¬ 
nel, i.e. the p-meson state. This choice uniquely fixes the 
pion mass dependence of the lepton mass to^(tops) and 
is subsequently used for all single-flavour contributions 
to the vacuum polarization function. 

H = i?phys = 1 reproduces the standard definition 
in Eq. 0 . Up to lattice artefacts, the standard defini¬ 
tion is also recovered at the physical value of the pion 
mass when the ratio 7J/7?phys becomes one 


Here u stands for the up quark, d for the down quark, 
c denotes the charm quark, and s the strange quark. 
Eq. ([^ shows that n^^{Q) results from the Fourier 
transformation of the correlator of two such currents. 
Taking up and down quarks together, since they are 
mass-degenerate in our setup, we decompose the quark- 
connected part of the hadronic vacuum polarisation 
tensor according to 

n,.AQ) = n;t{Q) + + i7;,(Q). (4) 


lim (tops) = (to^) . 

Henceforth we always use the definition in Eq. ([^ with 
H = my and drop the bar on the label for the lepton. 

The weight function w is known from QED pertur¬ 
bation theory as 


r2(l -I- ^Jl + A/rY^l + 4/r 


In our lattice calculation this decomposition into 
flavour contributions is particularly straightforward, be¬ 
cause for all quark flavours we use the one-point-split 
vector currents, which are conserved at non-zero lattice 
spacing and thus do not require further multiplicative 
or additive normalisation. From eq. Q, we can step¬ 
wise add the flavour contributions which will be done 
in the sections below. 

The standard integral definition in Eq. Q results in 
a strong non-linear pion mass dependence, in particular 
for the light quark contribution to ■ This behaviour 
originates from the introduction of the lepton mass to/ 
as an external scale, which is not related to the lattice 
parameters and in particular does not have an inherent 
value in lattice units. Employing Eq. 0 in the lattice 
calculation requires the input of the dimensionless com¬ 
bination a ■ mi and this renders the initially dimension¬ 
less quantity effectively dependent on the lattice 
scale setting. In view of this external scale problem. 


It has a pronounced peak at rpeak = Opeak/*^? = VS — 
2. As an illustrating example the corresponding peak 
locations for the electron, muon and tau are shown by 
the labels on the upper a;-axis in the upper plot of Fig. 
[T]for ensemble D30.48 in Table [T] below. 

For a thorough description of the lattice calculation 
and a proof of automatic 0{a) improvement of the vac¬ 
uum polarisation function we refer to [15j and |26j , re¬ 
spectively. In order to discuss systematic uncertainties 
later on, we briefly summarise our method of fitting the 
hadronic vacuum polarisation function here. 

First, the lowest lying vector meson masses, to/, and 
decay constants, ft, are determined from the time de¬ 
pendence of the two-point function of the light, strange, 
and charm point-split vector current, individually, at 
zero spatial momentum. Then n{Q‘^) determined in the 
momentum range between 0 and Q^ax is split into a 
low-momentum part for 0 < < Q^atch ^ high- 

momentum one for Q^atch < ^ Qmax is fitted 









4 


separately for each flavour and each ensemble. The low- 
momentum fit function is given by 

7t,ow(Q^) = E + E > (7) 

2=1 ^ j —0 

and the high-momentum piece is parametrised as fol- 
lows 

B-1 C-1 

= log(Q2) E + E • (8) 

fc=0 /=o 

They are combined according to 

7T(Q2) = (1 _ e(Q2 _ 

-\-0{Q — ) ) ( 9 ) 

where 0(x) is the Heaviside function. 

This defines our so-called MNBC fit function. Our 
standard fit for the light and strange quark contribu¬ 
tions is M1N2B4C1 which means M = 1, N = 2, 
i? = 4, and C = 1 in Eqs. Q and Q above. As 
value of Qmatch Heaviside functions in Eq. § 

we have chosen 2 GeV^. We have checked that varying 
the value of between f GeV^ and 3 GeV^ does 

not lead to observable differences as long as the tran¬ 
sition between the low- and the high-momentum part 
of the fit is smooth. For the upper integration limit we 
use Qmax = 100 GeV^, since the integrals are saturated 
there as can be seen in Fig. below. 

For each ensemble and flavour we perform a fit for 
n as given in Eq. From this ht we obtain the 
corresponding 17(0) and thus the subtracted polarisa¬ 
tion function. The latter is integrated using Eq. 
and the contributions from individual quark flavours 
are summed including the appropriate charge factors 
-I- e^, Cg, and e^. This results in an estimate for 
the hadronic leading-order lepton anomalous magnetic 
moment for each gauge field ensemble, which depends 
on the lattice spacing, the pion mass, and the lattice 
size, (a, mp 3 , L). As final step we per¬ 

form a combined extrapolation to the continuum and 
to physical quark masses. For this extrapolation some 
dependencies turn out to be negligible. The strange and 
charm quark mass in the sea and valence quark action 
have been tuned for each ensemble, so we do not need 
to consider these dependencies explicitly in More¬ 
over, the detailed discussion of the lattice data in the 
following sections will show that we only find significant 
lattice artefacts in the strange and charm quark contri¬ 
bution to and we will use an appropriately adapted 
fit ansatz. The dependence on the finite volume is dis¬ 
cussed in detail as part of the systematic error analyses 
in sections |3.3.1| and |4.3.1| and thus not part of the ex¬ 
trapolation described above. 


In [T71I28] the usage of Bade approximants has been 
advocated for fits in the small momentum region. The 
Fade fit functions are formally identical to the MN 
series of fits. We analysed the Fade approximants for 
our data in [5S]. We found agreement for the location 
of the poles provided the same number of fit parame¬ 
ters were used in both cases (See also the more elabo¬ 
rate study for the case of the muon performed by the 
RBG-UKQGD collaboration in [SO]-). We can thus em¬ 
ploy the same procedure as for the muon and show that 
it produces results compatible with phenomenological 
determinations for both the electron and the r-lepton 
without any modification. 

Our analysis has been performed on the same set 
of gauge field configurations [HES] as have been used 
in our previous work m- A detailed list of the lattice 
parameters can be found in table below. 


Ensemble 

a[fm] 

mps [MeV] 

L[fm] 

mps ■ L 



D15.48 

0.061 

227 

2.9 

3.3 

D30.48 

0.061 

318 

2.9 

4.7 

D45.32sc 

0.061 

387 

1.9 

3.8 

B25.32t 

0.078 

274 

2.5 

3.5 

B35.32 

0.078 

319 

2.5 

4.0 

B35.48 

0.078 

314 

3.7 

5.9 

B55.32 

0.078 

393 

2.5 

5.0 

B75.32 

0.078 

456 

2.5 

5.8 

B85.24 

0.078 

491 

1.9 

4.7 

A30.32 

0.086 

283 

2.8 

4.0 

A40.32 

0.086 

323 

2.8 

4.6 

A50.32 

0.086 

361 

2.8 

5.1 

cA2.09.48 

0.091 

135 

4.4 

3.0 


Table 1 Parameters of the Nf = 2-|-l-|-l flavour gauge 
field configurations that have been analysed in this work, a 
denotes the lattice spacing (cf. 1311 1. mps the value of the 
light pseudoscalar meson mass (cf. 1181 1 and L the spatial 
extent of the lattices. The right-most column gives the value 
for mps ■ L. The ensemble in the last line has Nf = 2 and 
physical pion mass. It is described in Refs. I20II21II32I . The 
ensemble name in the first column gives the bare quark mass 
in lattice units as the first pair of digits times 10“^ and the 
spatial lattice size L/a as the second pair of digits. 


3 The electron {g — 2) 

The LO hadronic contribution to the electron anoma¬ 
lous magnetic moment Ce is dominated by momenta 
below 10“^ GeV^. To a good approximation it can even 
be determined from the slope of the vacuum polarisa¬ 
tion at zero momentum oc dn/dQ^{Q^ = 0). There¬ 
fore, we only use the low-momentum part, 77iow(Q^), 
of the hadronic vacuum polarisation function Eq. ([^. 
The saturation of the integral for one of our ensembles. 
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namely B55.32 featuring tops « 390 MeV, a « 0.08 fm 

and L = 2.5 fm, is shown in the upper plot of Fig. for 

all three leptons by plotting 

„hvp/^2 
Wm, 


^i(Omax) = 


J 


( 10 ) 


af''P(100GeV^) ’ 

where is the LO hadronic contribution to 

the lepton anomalous magnetic moment integrated up 
to Qmax- This plot also implies that for the electron we 



OLax [GeV^l 



Fig. 1 Upper plot: Comparison of the dependence on the 
upper integration bound in Eq. n of the four-flavour lepton 
anomalous magnetic moments. The blue curve represents the 
ratio defined in Eq. ( |10[ | for the electron, the orange one for 
the muon, and the dark red one for the tau. i denotes 

the momentum value where the kernel function in Eq. fit 
attains its maximum. Lower plot: Comparison of the depen¬ 
dence of the relative statistical uncertainties of the integrands 
in Eq. fit on the squared momenta scaled by the lepton 
masses. The plots are based on data for the D30.48 ensemble 
featuring a = 0.061 fm, mps = 318 MeV, and L = 2.9 fm. 


have to rely mostly on the extrapolation of our vacuum 
polarisation data to the small momentum region. Al¬ 
though saturating well beyond momenta of O (l GeV^), 
also for the t, the renormalised vacuum polarisation 
function requires the subtraction of F?(0), which is de¬ 
termined from the same extrapolation to the small mo¬ 
mentum region. 


Despite the different masses of the leptons, the lower 
plot of Fig. shows that the relative statistical uncer¬ 
tainties are of the same orders of magnitude for all three 
leptons and display a universal dependence on /mf. 


3.1 Gontribution from up and down quarks 


The light quark contribution is depicted in Fig.[^ Here, 
we compare with the result at the physical value of 
the pion mass obtained with the standard definition 
Eq. Q on one ensemble iiniEi] with only one lattice 
spacing. 

For the extrapolation in the upper plot of Fig. we 
initially use 0 “*^ from all our ensembles, i.e. all volumes, 
lattice spacings, and pion masses. As the figure shows, 
with the present accuracy of the data we do not resolve 
any statistically significant dependence of on the 
lattice spacing or the lattice size. Using the modified 
definition (Eq. with H = my the ansatz 

(a, TOps, a) = A-b H TOpg (11) 


is then already sufficient to describe our data. The re¬ 
sult of this extrapolation is shown as the light grey band 
in Eig.[^ In principle, we are to add terms of higher or¬ 
der in TOpg for the extrapolation formula in Eq. (11), 


(a, rripg, L) = A + B TOpg + B 2 TOpg -(-... (12) 

to account for any non-linear dependence on rripg. The 
dark grey band in the background in Fig. shows the 
extrapolation with an additional term B 2 rripg. The dif¬ 
ference between this extrapolation and the previous one 
linear in the squared pion mass is insignificant. This in¬ 
significance of terms in TOpg beyond the linear one when 
using the improved definition in Eq. (§ turns out to be 
a universal property of a( for all leptons and all flavour 
combinations / = ud, udsc considered here and below. 
In addition, we check explicitly for lattice artefacts in 
the extrapolation by adding a term C a? to Eq. ( |lT| ) 
with the result shown in the lower plot of Eig. 

We assume that lattice artefacts for the data at the 
physical point are negligible as well, which gives merit 
to the observed compatibility of the physical point re¬ 
sult with the value determined by the linear extrapola¬ 
tion of the data described here. 


3.2 Adding the strange and the charm quark 
contributions 

When incorporating the heavy, second-generation flavours, 
which are described by the Osterwalder-Seiler action [321 
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mis [GeV^ 



Fig. 2 Upper plot: light-quark contribution to with 

filled symbols representing points obtained with Eq. ([^, open 
symbols refer to those obtained with Eq. 0 , i. e. H = 1. In 
particular, the two-flavour result at the physical point has 
been computed with the standard definition. The light grey 
errorband belongs to the linear fit, whereas the dark grey 
errorband is attached to the quadratic fit. Lower plot: com¬ 
bined chiral and continuum extrapolation of the light-quark 
contribution to 0 ^*° allowing for lattice artefacts. 


134) and whose masses have been tuned to their physical 
values as shown in [TS], we take 0{a^) lattice artefacts 
into account. The four-flavour result for at the 

physical point in the continuum limit is obtained from 
simultaneously extrapolating in the pion mass, mps, 
and to zero lattice spacing a using the ansatz 

(mps, a) = A-k i?TOpg -I- Ca^ . (13) 

A, B, C denote the free parameters of the fit. In the 
presence of the strange and charm contributions to 
the parameter C will also contain terms ^ Wg p, p 
from the renormalised charm and strange quark mass 
and also receives contributions from lattice artefacts 
possibly present in the light quark contribution to Og'^P. 
The reason for omitting a linear term in a is that auto¬ 
matic 0{a) improvement is retained for our definition 
of the hadronic vacuum polarisation function at maxi¬ 


mal twist [26] . As we have discussed in [15] , systematic 
effects from varying the heavy valence and sea quark 
masses within the range given there have been found 
to be negligible. This is partly due to the contribution 
of strange and charm quark current correlators to vac¬ 
uum polarisation being at least an order of magnitude 
smaller than those from the light quarks. 

For Og'^P the corresponding fit is shown in Fig. 
Again we use data for ag''P for all lattice spacings, pion 
masses, and lattice volumes in this extrapolation. The 


ansatz in Eq. (13) is sufficient for a good description of 


all our data. This is shown by the three dashed lines, 
which evaluate the fit function for the individual lat¬ 
tice spacings. We have checked that amending the fit 
function by higher powers in Wpg does not lead to sig¬ 
nificantly different results for the extrapolated value. 
Comparing the results for different lattice volumes for 
lattice spacing a = 0.078fm in Fig. [^suggests the ab¬ 
sence of observable finite volume effects. However, for 
the compilation of our complete error budget we inves¬ 
tigate these effects in more detail below. 

Our result with only statistical uncertainty is 


z^'^P = 1.78 (06) • 10 


-12 


(14) 


2.0e-12 


1.8e-12 



1.6e-12 





1.4e-12 - 
1.2e-12 - 

l.Oe-12 - 
8.0e-13 - 

6.0e-13 -^^— 

0 ml 0.05 


a ^ 0 result ■ t ^ 

dispersive analysis «—B—« 
a = 0.061 fin, L = 1.9 fm >—• 
a = 0.061 fin, L = 2.9 fm >—■—> 
a = 0.078 fm, L = 1.9 fm « — a —i 
a = 0.078 fm, L = 2.5 fm «—•—« 
a = 0.078 fm, L = 3.7 fm «— ■ —« 
a = 0.086 fm, L = 2.8 fm •—♦—> 

aefmps, 0.061 fm) . 

aefmps, 0.078 fm) . 

aefmpsi 0.086 fm) --- — 


0.1 0.15 0.2 0.25 

mps [GeV^] 


Fig. 3 Chiral and continuum extrapolation of the Nf = 
2 -k 1 -k 1 contribution to . The inverted red triangle shows 
the value extrapolated to the continuum and to the physical 
value of the pion mass. It has been displaced to the left to fa¬ 
cilitate the comparison with the dispersive result in the black 
square [35]. 


3.3 Systematic uncertainties 


In this section, we give an account of systematic un¬ 
certainties of our result for Og'^P given in Eq. (14). We 
have investigated finite size effects (FSE), the depen¬ 
dence of our chiral extrapolation on the incorporation 
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of large pion masses, vector meson fit ranges, and the 
dependence of our results on different vacuum polarisa¬ 
tion fit functions. Moreover, for one ensemble the light 
quark-disconnected contribution is quantified. 


3.3.1 Finite size effects 

As described in detail in Ref. m, the Nf = 2 -1-1 -I-1 en¬ 
sembles analysed in this work feature 3.35 < mps L < 
5.93, where L is the spatial extent of the lattice. Re¬ 
stricting our data to the condition mpsL > 3.8 and 
Tops A >4.5, respectively, yields 

[mpsL > 3.8) = 1.77 (07) • 10"^^ , (15) 

(mpsi > 4.5) = 1.83 (10) • lO'^^ , (16) 


after combined chiral and continuum extrapolation. This 
matches the result given in Eq. (141 and thus indicates 
that FSE are negligible in our computation. This find¬ 
ing is supported by comparing the results of two en¬ 
sembles only differing in lattice size provided in Tab. 
The numbers do not change when restricting the mo¬ 
menta of the larger ensemble to those of the smaller 
one. The FSE attributed to the lowest achievable mo¬ 
mentum being ^ mixes with FSE entering the choice of 
different fit functions. We take a conservative approach 
and consider these effects separately. 


Ensemble (^)^xT 

\ a / a 

hvp hvp 

“e,ud “e P 


B35.32 323 ^ 64 

B35.48 483 x 96 

1.44(05) • 10-12 1.66(05) • IO -12 

1.44(05) • 10-12 1.69(05) • IO -12 


Table 2 Comparison of light-quark contribution to 
and total Oe'^P from ensembles of different volumes. 


3.3.2 Chiral extrapolation 

We have checked the validity of the chiral extrapola¬ 
tion by restricting the data, comprising pion masses be¬ 
tween 227 MeV and 491 MeV, to the condition mps < 
400 MeV. The value we obtain 


3.3.3 Vector meson fit ranges 

Our standard computation involves the determination 
of the masses and decay constants of the vector me¬ 
son ground states for the different flavours. Their val¬ 
ues depend on the choice of fit ranges. We have anal¬ 
ysed different ht ranges for the two-point functions of 
the light, strange, and charm vector currents and prop¬ 
agated the uncertainties to the values for oIC^p. This 
showed that excited state contaminations are significant 
only for my and fy determined from the light vector 
current-current correlator. Variations of the standard 
fit ranges by 0.1 fm to the left, right and both simul¬ 
taneously do not lead to any observable differences in 
Og'^P for the sy^s- and the J/f) correlator. Furthermore, 
the heavy flavour contributions are approximately one 
order of magnitude smaller than the light quark contri¬ 
bution such that their systematic uncertainties would 
not noticeably impact the overall uncertainty of a^'^P. 

In the upper plot of Fig. the dependence of the 
light quark contribution to the electron anomalous mag¬ 
netic moment on the fitrange for the p-correlator is plot¬ 
ted. The lower limit 0.6 fm is a lower bound for the time 
region, where a single-state fit describes the 2-point cor¬ 
relation function. The upper end at 1.3 fm is less strin¬ 
gent, but the signal to noise ratio for the masses dete¬ 
riorates quickly and the differences in a plot like Fig. 
become insignificant. 

Taking half the difference of the central values ob¬ 
tained for [0.6 fm, 1.2 fm] and [0.7 fm, 1.2 fm] gives a sys¬ 
tematic uncertainty of 

Ay = 0.035 • 10-1^ . (18) 


3.3.4 Number of terms in MN fit function 

The number of terms in the fit function Eq. Q is given 
by M and N. M1N2 is our standard choice. Repeating 
the whole analysis with different numbers of terms for 
the light quark contribution leads to the results shown 
in the lower plot of Fig. [^ We observe that the chirally 
and continuum extrapolated results of fit functions in¬ 
volving one and two poles are not compatible and thus 
we assign a systematic error by taking half the differ¬ 
ence of the central values of the result of the M2N3 and 
the M1N2 fit. This leads to a systematic uncertainty of 


x[)''P = 1.78 (07) • 10 


-12 


(17) 


only features a slightly larger uncertainty compared to 
the result in Eq. (14). Thus, we do not assign a sys¬ 
tematic uncertainty to the usage of pion masses above 
400 MeV. 


= 0.071 


10 


-12 


(19) 


This results in the dominant systematic uncertainty of 
the determination of Og'^P. For the strange quark the 
systematic uncertainty from different values of M and 
N is 


= 0.007 • 10 


-12 


( 20 ) 
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'-B-' 


[0.7,1.3]fin 

'-B-' 


[0.7,1.2]fiu 


1—B—1 

[O.e, 1.3]fni 


1—B—1 

[O.e, 1.2]tiu 


>12 

1.5e-12 

1.6t>12 

1.7{^12 



'-B-' 

M2N4 



■ □ ' 

M2N3 



'-B-' 

M2N2 


>—B—' 


M2N1 


'-B- 


M1N4 


'—B—' 


M1N3 


' □ ' 


M1N2 


1.4f>12 1.5e-12 1.6f>12 1.7e-12 1.8e-12 


Fig. 4 Dependence of on the fitrange of the p-correlator 
(upper plot) and on values chosen for M, N in the vacuum po¬ 
larisation fit function (lower plot). The standard p-correlator 
fit range is [0.7 fm, 1.2 fm] and the standard fit function cor¬ 
responds to M1N2. 

which we add to the light quark one. The differences of 
results from different fit functions for the charm quark 
contribution have turned out to be negligible such that 
the total systematic error originating from employing 
various numbers of terms in the fit function amounts to 

AmN = 0.078-10-^^ . ( 21 ) 

3.3.5 Disconnected contributions 

Leaving out the quark-disconnected contributions is a 
systematic uncertainty we cannot completely quantify, 
yet. We have started investigating their magnitude on 
the B55.32 ensemble mentioned already before. Using 
the local vector current we have detected a signal for 
the light quark part of the vacuum polarisation function 
when using 24 stochastic volume sources on 1548 con¬ 
figurations and 48 stochastic volume sources on 4996 
configurations. Employing the one-end trick [37], the 
isovector part 

n%{x,y) = {Jl{x)Jl{y)) ( 22 ) 

with is significantly different from zero. 

However, this is a pure lattice artefact and will not 



Fig. 5 Comparison of the light quark contributions to the 
unsubtracted hadronic vacuum polarisation function from 
quark-connected and disconnected diagrams of the local cur¬ 
rent correlator. Zy has been obtained from the ratio of the 
connected part of the conserved and local current-current cor¬ 
relators. The values have been computed with the anal 3 rtical 
continuation method described in |.36| without correcting for 
finite-size effects. 

contribute in the continuum limit. On the other hand, 
the more interesting isoscalar part 

nl^{x,y) = ^(J°(a;)J°(y)) (23) 

with = \xiu^X is compatible with zero. The con¬ 
nected and disconnected pieces of the polarisation func¬ 
tion for the light flavours are depicted in Fig. 

A comparison of the values of for all three lep¬ 
tons on the B55.32 ensemble with and without incor¬ 
porating the disconnected contributions is presented in 
Tab. Here, we have combined the connected pieces 
obtained from the point-split current correlator with 
the isoscalar part of the disconnected contributions ob¬ 
tained from the local current correlator using the renor¬ 
malisation constant Zy determined from the ratio of 
the connected pieces of the conserved and the local 
vector current two-point functions. Therefore and be¬ 
cause we only have results for one ensemble, the num¬ 
bers below can only give hints on the influence of the 
disconnected pieces. We observe the tendency that for 
all three leptons decreases when incorporating the 
disconnected contributions as has been predicted in [ 35 ] ■ 
However, this is statistically not significant. Further¬ 
more, we find that the magnitude of the disconnected 
contributions is comparable to our current uncertainty. 
Hence, it will be mandatory to compute them when 
aiming at more precise results. For the muon the value 
shifts by « 3%, which is also not statistically signifi¬ 
cant at this stage, but is in accordance with the upper 
bound of 4 — 5% given in [ 33 ] as well as more recent 
high-statistics evaluations in iiniiii]. 
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I without disc with disc 


“Ud 

1.44(04) • 10-12 

1.39(07) • 10-12 

hvp 

“/j,ud 

5.42(14) • 10"® 

5.26(25) • 10-® 

hvp 

“r.ud 

1.27(03) • 10-® 

1.24(04) • 10-® 


Table 3 Comparison of light-quark contributions to 
with and without disconnected pieces in the low-momentum 
region for the B55.32 ensemble. For all contributions the re¬ 
definition Eq. and our standard analysis have been used. 


The disconnected heavy flavour contributions need 
to be considered as well. We plan to check their size 
in future calculations. The pure charm quark contribu¬ 
tions have been computed in perturbation theory and 

shown to be suppressed by a factor [l 2 ], where 

is the relevant energy scale of the problem. 

3.4 Comparison with the phenomenological value 

Adding the quantified systematic uncertainties in quadra¬ 
ture we obtain as final result 

= 1.782 (64)(86) • ^ ( 24 ) 

This can directly be compared with the phenomenolog¬ 
ical determination of |35] 

a^’^P = 1.866 (10) (05) • lO’^^ ^ (25) 

They are fully compatible with each other although our 
lattice result still is afflicted with larger errors. 


4 The T-lepton {g — 2) 

The large mass of the tau lepton, nir ~ 1.8 GeV, im¬ 
plies a peak of the weight function in the expression 
for the LO hadronic contribution to its magnetic mo¬ 
ment in Eq. Q at = 0.745 GeV^. This is very 

different from the peak position of the electron weight 
function. Hence, the saturation of a(l''P requires data 
from a different part of the subtracted vacuum polari¬ 
sation function, in particular, also the high-momentum 
piece of our fit function Eq. (|^ is important here. 


4.1 Contribution from up an down quarks 

As for the electron, we start off by showing the contribu¬ 
tion of the first-generation flavours to a(l''P in the upper 
plot of Fig.[^ The data show a qualitatively similar be¬ 
haviour to those of the electron in Fig. Their values 


differ, however, by six orders of magnitude. In partic¬ 
ular, by comparing the upper and lower plot of Fig. 
we And that no significant lattice artefacts are present 
and that the data at unphysical pion masses obtained 
with Eq. (§, can be linearly extrapolated to the physi¬ 
cal point. This demonstrates again that the method of 
including in the weight function is advantageous 

^“phys 

for the chiral extrapolation. The value extrapolated in 
this way and using all available lattice ensembles agrees 
with our calculation directly at the physical pion mass 
shown as the open square in Fig. 



2.80e-06 

2.60e-06 

2.40e-06 

2.20e-06 

2.00e-06 

1.80e-06 

1.60e-06 

1.40e-06 


0 ml 0.05 0.1 0.15 0.2 0.25 




I 4 Jt 


a —0 result 
Nf = 2 result 
a = 0.061 fni, L = 1.9 fin 
a = 0.061 fm, L = 2.9 fm 
a = 0.078 fni, L = 1.9 fin 
a = 0.078 fni, L = 2.5 fin 
a = 0.078 fm, L = 3.7 fm 
a = 0.086 fm, L = 2.8 fin 
a^'^(mps, 0.061 fm) 
a^'^Jmps,0.078 fm) 
(mps, 0.086 fm) 




mis [GeV^] 


Fig. 6 Upper plot: light-quark contribution to with 

filled symbols representing points obtained with Eq. , open 
symbols refer to those obtained with Eq. 0 , i. e. H = 1. We 
note that the two-fiavour result at the physical point has been 
computed with the standard definition. The light grey error- 
band belongs to the linear fit (dotted black line), whereas the 
dark grey errorband is attached to the quadratic fit (solid 
green line). Lower plot: combined chiral and continuum ex¬ 
trapolation taking into account leading order lattice artefacts. 
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4.2 Adding the strange and the charm quark 
contributions 


As for the electron, we perform the chiral and contin¬ 
uum extrapolation of the complete four-flavour result 
using a fit of the form given in Eq. and data for 
all lattice spacings, pion masses, and lattice volumes 
simultaneously. It is shown in Fig. Comparing this 
with Fig. 1^ we see that the lattice artefacts are much 
smaller than for the electron such that we would have 
obtained a compatible result when omitting the term 
in Eq. (13). As can be seen in Figs. |^and|9, for the tau 
lepton both, the strange and the charm contribution 
do not show significant cut-off effects and hence, also 
for the total contribution effects are small. We never¬ 
theless perform the continuum extrapolation in order to 
use exactly the same analysis strategy as for the other 
leptons. 


3.5e-06 


3.0e-06 

I 

3 b- 
Oj 

2.5e-06 


2.0e-06 


0 ml 0.05 0.1 0.15 0.2 0.25 

mps [GeV^] 

Fig. 7 Chiral and continuum extrapolation of the Nf = 
2 -|- 1 -h 1 contribution to The inverted red triangle 

shows the value in the continuum limit at the physical value 
of the pion mass. It has been displaced to the left to facilitate 
the comparison with the dispersive result depicted as black 
square 1431 . 




a —> U result •— 

dispersive analysis •—B—• 
a = 0.061 fnn L = 1.9 fm — 
a = 0.061 fm, L = 2.9 fm •—■—• 
a = 0.078fm, L = 1.9fm — 
a = 0.078fm, L = 2.5fm i 9 ■ 
a = 0.078fm, L = 3.7fm 1—■—1 
a = 0.086 fm, L = 2.8 fm 1 — «—1 

aefmpg,0.061 fm) . 

aefmps,0.078fm) . 

ae(7nps'0.086fm) . 


Our four-flavour result with only statistical uncer¬ 
tainty reads 


= 3.41 (8) • 10-® . 


(26) 


4.3 Systematic uncertainties 

We have investigated the same systematic uncertainties 
for our determination of a(l''P as for the case of the elec¬ 
tron. The influence of the disconnected contributions 
has already been discussed in the section of 



Fig. 8 Continuum limit of strange quark contribution to 
at approximately fixed pion mass. 

5.0e-07 -,- 


4.5e-07 


U 4.0e-07 


3.5e-07 - 0 


3.0e-07 


data at lixed mpg 320 MeV^ 
CL with linear fit 


0.002 0.004 0.006 0.008 0.01 


Fig. 9 Continuum limit of charm quark contribution to 
at approximately fixed pion mass. 


4-3.1 Finite size effeets 


Restricting our data to the conditions mpsA > 3.8 and 
mpgL > 4.5 yields 


a^'^P (ttipsA > 3.8) = 3.40 (09) • 10"® , (27) 

a^'^P {mpsL > 4.5) = 3.54 (13) • 10"® . (28) 


This is compatible with the result in Eq. (26). Compar¬ 
ing again the two ensembles at mpg « 315 MeV which 
only differ in the extent of the lattices also indicates 
negligible finite size effects as shown in Tab. Hence, 
we do not assign a FSE related systematic uncertainty. 


Ensemble ( — ) x — 

V a / a 

hvp hvp 

B35.32 323 ^ 64 

B35.48 483 x 96 

2.62(06) • 10-6 3.40(07) • 10-6 

2.60(067 10-6 3.41 (07) • 10-6 


Table 4 Comparison of light-quark contribution to a^''^ 
and total from ensembles of different volumes. 
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4-3.2 Chiral extrapolation [0.7fm, 1.2 fm] results in an estimated systematic un¬ 

certainty of 

Restricting the analysed ensembles to those featuring 

pion masses mps < 400 MeV, we get 4\y = 0.046 • 10 ® . (30) 


= 3.45 (09) • 10"® . 


(29) 4-3.4 Number of terms in MNBC fit function 


This is again compatible with the value given in Eq. (26). 
Hence, we do not assign a systematic uncertainty to the 
fact that ensembles with pion masses above 400 MeV 
have been employed when extrapolating to the physical 
value of the pion mass. 


4-3.3 Vector meson fit ranges 

The situation is similar to the case of the electron re¬ 
ported above. Only the excited state contamination in 
the p-correlator has to be taken into account as sys¬ 
tematic uncertainty. In the upper plot of Fig. the 
dependence of the light quark contribution, on the 
fit range chosen to extract the spectral information from 
the p-correlator is depicted. 


Due to the large Qpgajj we have to take the whole vac¬ 
uum polarisation function Eq. into account, includ¬ 
ing in particular the high-momentum piece in Eq. Q. 
Thus, we have four different types of terms in the fit 
function that can have different numbers of summands. 
We only find observable differences in the light quark 
sector. But also here the results from different fits are 
all compatible as shown in the lower plot of Fig. [T^ 
Conservatively, we take half the difference between the 
M2N3B4C1 and M1N2B4C1 ht and assign a systematic 
uncertainty of 

Amnbc = 0.032 • 10-6 (31) 

to our choice of the fit function. 


1 —a— 1 


[0.7,1.3]fm 

' -B-1 


[0.7,1.2]fm 


1 —a— 1 

[0.6,1.3]fm 


1 —a— 1 

[0.6,1.2]fm 


2.5fr06 2.6fr06 2.7fr06 2.8^06 2.9e-06 3.0c-06 


'-B-' 

M2N4B4C1 

'-B-' 

M2N3B4C1 

'-B-' 

M2N2B4C1 

'-B-' 

M1N4B4C1 

'-B-' 

M1N3B4C1 

'-B-' 

M1N2B4C2 

'-B-' 

M1N2B4C1 

'-B-' 

M1N2B3C2 

'-B-' 

M1N2B3C1 


2 . 60 O 6 2.7o06 2.8e-06 2.9o06 3.0o-06 


Fig. 10 Dependence of a'^'^ on the fit range of the p- 
correlator (upper plot) and on the values chosen for M, N, 
B, and C In the vacuum polarisation fit function (lower plot). 
The standard p-correlator fit range is [0.7 fm, 1.2 fm] and the 
standard fit function corresponds to M1N2B4C1. 


Taking half the difference of the central values ob¬ 
tained for [0.6 fm, 1.2 fm] and our standard fit range 


4.4 Comparison with the phenomenological value 

Including the identified systematic uncertainties added 
in quadrature, our final four-flavour result reads 

a6''P = 3.41 (8)(6 )-10-6. (32) 

This agrees with the one obtained by a dispersive anal¬ 
ysis [43] 

= 3.38 (4) • 10-6 . (33) 

Compared to the electron, even better agreement be¬ 
tween the lattice and the phenomenological result is 
observed for the r-lepton. In this case, the uncertainty 
of our twisted mass LQCD calculation is only about 
twice the phenomenological one. 


5 Summary and Conclusions 

In this article we have presented the first four-flavour 
LQCD computation of the LO hadronic vacuum polar¬ 
isation contributions to the anomalous magnetic mo¬ 
ments of the electron and the r-lepton. Our results 
have been obtained with Nf = 2 -|- 1 -|- 1 twisted mass 
fermions mostly at unphysically large pion masses but, 
at least for the light quark contribution, also directly at 
the physical point. We find that for both, the electron 
and the tau lepton, the chirally extrapolated values for 
the light quark contributions agree with the one at the 
physical point. 
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For our data at unphysically large values of the pion 
mass we have investigated the systematic uncertainties 
of the method used to obtain our final results. In par¬ 
ticular, we have addressed the effects of non-zero lattice 
spacings, the finite volumes, the ht range for extracting 
the vector meson mass, and using different fit functions 
for the vacuum polarisation function. As an additional 
uncertainty we have investigated the disconnected con¬ 
tributions on one of our 4-flavor ensembles (B55.32) by 
using the local vector current. This led to the first ob¬ 
servation of a signal for the disconnected diagrams dur¬ 
ing our calculations, which, however, is compatible with 
zero within our current errors and which we therefore 
have neglected. This will no longer be justified once the 
uncertainties of the connected pieces are reduced and a 
full quantification of the quark-disconnected contribu¬ 
tion will become significant. 

Our final results are summarised in Tab. below 
and agree with the phenomenological determinations of 
the electron and tau lepton magnetic moments which 
are also shown there. This universal agreement across 
all three leptons and thus distinct weightings of the 
subtracted polarisation functions is elucidated by our 
findings in Ref. [H]. There it was shown, that the sub¬ 
tracted vacuum polarisation function itself calculated 
with the methods used in this work and described in 
more detail in Ref. m, is compatible with the phe¬ 
nomenological result for in the range 0 < < 

C>(lOGeV2). 


I this work dispersive analyses 


hvp 

1.782 (64)(86) • IQ-^^ 

1.866(10) (05) • 10-12 


hvp 

6.78(24)(16) • 10"® 

6.91 (01) (05) • 10"® 

m 

hvp 

ar 

3.41 (8)(6) • 10-® 

3.38 (4) • 10-® 

m 


Table 5 Comparison of our first-principle values for ae'^^, 
ajl''®’, and with phenomenological results. 


As expected from the graph in the lower plot of 
Fig. a the relative statistical uncertainties in all three 
cases are similar. For the electron the systematic un¬ 
certainty already exceeds the statistical one. 

As in the case of the muon, also for the electron 
and tau lepton anomalous magnetic moments the errors 
of our calculations are still larger than those from the 
dispersive analyses quoted above. However, it can be 
expected that with futnre lattice QCD calculations at 
the physical value of the pion mass, increased statistics 
and an even better control over systematic uncertainties 
the phenomenological error can be matched, if not even 
beaten, especially for the r-lepton. 
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